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Recently,  several  authors  [1-3]  have  speculated  that  small  magnetic 
fluctuations  associated  with  drift  waves  may  be  responsible  for  anomalous 
electron  energy  losses  in  Tokamaks.  These  generally  extend  earlier  work 


[A]  on  the  effect  of  stochastic  wandering  of  field  lines  on  particle 
transport.  In  Refs.  [1-3]  the  basic  approach  is  to  examine  particle 
loss  rates  in  an  ambient  magnetic  field  with  a steady  random  [1-3]  or 
oscillating  [1]  radial  field  perturbation.  Ih  order  to  get  expressions 
for  diffusion  and  thermal  conduction,  these  authors  make  various  assump- 
tions about  magnetic  island  width  [1]  and/or  correlation  lengths  of  the 
field  fluctuations  [1-3]. 

Three  basic  conclusions  stem  from  this  work.  Firstly,  very  small 
magnetic  perturbations  can  have  a large  effect  on  transport.  Secondly, 
electrons  are  lost  along  the  field  lines  much  faster  than  the  ions,  due 
to  their  greater  velocity.  Thus  the  diffusion  is  not  naturally  ambi- 
polar  and  large  radial  electric  fields  must  be  set  up  to  confine  the 
3 electrons.  Finally,  particles  with  larger  parallel  velocity  are  lost 

much  faster  than  particles  with  smaller  parallel  velocity,  simply 
because  they  move  faster  along  the  field  lines.  Although  this  point 
was  not  emphasised  in  Refs.  [1-3],  it  is  clearly  pointed  out  in  related 
work  concerning  diffusion  by  randomly  wanderir  magnetic  fields  of 
energetic  cosmic  rays  in  astrophysical  plasmas  [ .*>]  , and  diffusion  of 
energetic  electrons  in  laser  produced  plasmas  [6]. 

While  the  first  result  is  certainly  very  interesting,  the  second  is 
somewhat  unsettling  and  the  third  is  not  borne  out  by  experiment.  Indeed, 
experiments  generally  show  that  the  confinement  time  of  runaway  electrons 
in  a Tokamak  plasma  is  comparable  to  or  longer  than  the  overall  energy 
confinement  time  [1]. 

The  purpose  of  this  note  is  to  show  that  ordinary  quasi-1 inoar 
Note:  Manutcript  submitted  February  22,  1978. 
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theory  [8,  9]  gives  a quantitative  description  of  particle  diffusion 
and  thermal  conduction  associated  with  magnetic  flutter  of  drift  waves. 
It  is  not  entirely  clear  whether  quasi-linear  theory  is  valid;  however 
it  is  not  clear  that  it  is  invalid  either.  Refs.  [1-3]  all  invoke 
some  sort  of  coherent  or  steady  random  structure  of  the  field.  For 
instance  Ref.  [1]  assumes  a single  drift  wave  on  a rational  surface  and 
invokes  an  assumed  island  width.  Refs.  [2  and  3]  assume  steady  magnetic 
field  fluctuations.  However,  if,  at  each  point,  the  turbulence  has  a 
wide  spectrum  of  wave  numbers  and  frequencies,  as  seems  to  be  indicated 
by  scattering  experiments  on  ATC  [10]  and  TFR  [11],  it  may  well  be  that 
ordinary  quasi-linear  theory  is  the  more  accurate  description. 

In  any  case,  quasi-linear  predictions  of  transport  induced  by 
magnetic  flutter  are  simply  derived  for  a known  spectrum.  This  alone 
argues  that  the  quasi-linear  expressions  should  be  recorded.  A quasi- 
linear  treatment  also  has  the  incidental  advantage  that  it  is  not 
necessary  to  assume  "cit,jdic"  wandering  of  the  magnetic  field  lines. 

We  show  here  that  the  quo-i-linear  expressions  do  indeed  indicate 
that  small  magnetic  flutter  can  enhance  the  electron  energy  flux. 

However  they  also  show  that  the  particle  flux  is  ambipolar  and  thac 
there  is  no  enhanced  transport  of  runaway  electrons,  which  do  not 
diffuse  - simply  because  they  cannot  be  resonant  with  any  wave  of  the 
spectrum. 

We  also  show  how  particle  and  energy  fluxes  arc  computed.  Our 
calculation  closely  parallels  that  given  in  Ref.  [8].  We  begin  by 
assuming  a steady  state  and  cylindrical  geometry,  but  with  arbitrary 
theta  and  z components  of  the  average  magnetic  field.  The  averaged 
electron  Vlasov  equation  is 


.... 


■ :&■[(«> + 7^ x «>)  <V  * ((*  * * *)  ?.)] 


“ 0 


(1) 

where  < > indicates  an  ensemble  average,  tildes  indicate  fluctuating 

parts,  f is  the  electron  distribution  function  and  the  other  notation 
e 

is  standard.  In  Equation  (1)  collisional  drag  between  electrons  and 
ions  is  neglected  as  it  does  not  contribute  directly  to  wave  induced 
particle  flux.  We  multiply  equation  (l)  by  v and  integrate  over  v 
to  obtain 


e <£>  <ne>+  e <E  ne>  + 


~ J X <fe>  x <B>  + ~ ^ J&'v  - X 


(2) 


where  n^  is  the  number  density  of  electrons.  In  equation  (2)  the 
contribution  from  the  vf  3<fe>/0r  term  has  been  ignored  as  it  can 
easily  be  shown  [8]  to  be  smaller  than  the  contribution  from  the 
(v  X <B>)<fe>  term  by  the  ratio  of  the  electron  Larmor  radius  to  the 
macroscopic  scale  length. 

We  now  take  the  0-component  of  equation  (21.  As  the  plasma  is  in 
a steady  state?  <E^>  = 0 , as  indicated  by  Maxwell's  homogeneous  curl 
equation.  The  remaining  terms  are 

e <W  - c <V  *>fl  = ! Bo/div  vr  <fe>  * (3) 

where  = B^z  and  J_e  , the  fluctuating  electron  current  density 

is  defined  by 

s - ey*dJv  v . (A) 

The  integral  on  the  right  hand  side  of  equation  (3)  is  just  the 
radial  electron  particle  flux,  , say.  Rearranging  equation  (3)  we 
find 
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An  analogous  expression  for  the  ion  flux  may  be  derived  by  the 
same  method  (provided  the  ratio  of  the  ion  Lartnor  radius  tc  the  macro- 
scopic scale  length  is  small).  It  is 


ri  ' * {<;i  T>  + eT  <h  * !>»]  • 


(6) 


Note  that  equations  (5)  and  (6)  are  the  standard  expressions  for  the 
particle  fluxes  [9],  augmented  by  a contribution  from  the  magnetic 
flutter  [1].  Since  the  plasma  is  quasi-neutral  (<n^>  = <n^>  anc* 
n^  = n^  ) it  is  obvious  that  the  particle  flux  arising  from  the  electric 
field  is  ambipolar.  It  is  not  obvious  that  the  particle  flux  arising 
from  the  fluctuating  magnetic  field  is  ambipolar,  since  in 

general.  However,  if  we  take  the  cross-product  of  Maxwell's  inhomogen- 
eous curl  equation  with  B and  average,  w<5  obtain 


Uir  ~ 


<B  x ( V x B)  > = — <B  X (Je  + J£)  > . 


(7) 


The  left  hand  side  can  be  reduced  to  the  divergence  of  the  electro- 
magnetic momentum  stress  tensor.  Since  there  is  no  variation  of  any 
averaged  quantity  in  the  6 or  z directions,  the  0-componer.t  of 
equation  (7)  reduces  to 


4~  r*  <3  B-  > = <B  x (J  +J.  )v.. 
r 2 dr  r 6 N-  ■ -e  — i 6 


Whatever  the  radial  eigenfunction,  it  does  not  seem  reasonable  to  us 
that  the  average  force  the  plasma  exerts  on  itself  can  be  balanced  by 
•—  r 1 <5^  . Therefore  x J ■>  = 0 and  the  particle  diffusion 
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Is  ambipolar. 


Thus  the  wave  frequency  and  polarisation  are  picked  out  so  that 

n = n.  and  (B  X J > = - <B  X J, \ . This  is  analogous  to  the  case  of 
e i - — e — — i 

purely  electrostatic  waves,  where  the  polarisation  is  known  but  where 
the  frequency  is  picked  out  so  that  ng  = n^  . 

We  now  turn  to  the  problem  of  the  anomalous  electron  energy  flux 
driven  by  the  magnetic  flutter  associated  with  drift  waves.  The  calcula- 
tion again  parallels  that  given  in  Ref.  [8].  Taking  the  imv:v,  moment 

o 

of  the  Vlasov  equation,  averaging  over  wave  phases,  and  neglecting  the 

[tt—  + v.V]  <f  > terms  on  the  same  grounds  as  in  our  treatment  of  the 
at  — — e 

particle  fluxes,  we  find 

/•  mv 2 v , f.  2 £ » ..  8f  \ , 9 <f  >1 

„ . i/d>.  —i  -57"j  <8> 


By  various  partial  Integrations,  the  last  term  in  equation  (81  can  be 
manipulated  into  eB/mc  times  the  radial  electron,  energy  flux,  Wg  . 

The  other  terms  in  equation  (8)  can  also  be  simplified  by  partial  inte- 
grations, so  that  equation  (8)  reduces  to 
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f > 
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<m  vev  . E 
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V 
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(9) 
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where,  in  deriving  equation  (9),  we  have  neglected  terms  like  — B^ 
because  we  are  only  interested  in  electron  energy  flux  driven  by  radial 
magnetic  flutter.  Thus  equation  (9)  is  an  expression  for  the  energy 
flux  in  terms  of  a given  wave  spectrum. 


Let  us  assume  that  linear  theory  gives  the  result 


fgCk,  vj  u)  = E(k)  . F(k,  Vj  ui)  (10) 

in  the  absence  of  magnetic  flutter.  In  the  presence  of  radial  magnetic 
flutter,.  E in  the  Vlasov  equation  is  replaced  by  E + ^ v x B , so 
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tlijt  in  the  presence  of  magnetic  flutter  : 


fe(k,  v,  u)  = ( E(k)  + ^ v x B(k))  . F(k,  v,  u>)  (11) 

Assuming  that  tlve  0-component  of  the  force  is  dominant  with  or  without 
magnetic  flutter  and  also  that  B^  = 0 , the  equation  for  energy  flux 

reduces  to 

we  = + ~C  Sr(^l )+mvKV^EV^4  Br(k))>J  F(k,  v, 


+ c.c. 


Obvious  generalisations  of  equation  (12),  which  relax  our  simplifying 
assumptions,  can  be  made. 

Kquation  (12)  show'<  the  energy  flux,  with  or  without  the  radial 
magnetic  flutter.  Clearly,  a measure  of  the  importance  of  magnetic 
flutter  is  simply  the  ratio  B^(k)/c  Eg(ki  . Adopting  the  result 
of  Ref.  [ 1]  : 

®r^-*  iu  fie  AjmM  ( Ti  U*e\  e®(k)  ,,, 

T"  iT-Vi-v  tt-HP"-  (l  + T13“J  • -r—  ' (13 

II  \ e / e 


one  can  show  that  for  Tokamak  parameters  v B/cF.  can  Indeed  be 
of  order  unity  or  larger. 


As  is  apparent  from  equation  (12),  only  the  real  part  of 

F(k,  v,  u)  contributes  to  the  energy  flux.  Depending  on  whether  a 

particle  is  magnetically  trapped  or  not,  the  real  part  of  F(k,  v,  j) 

goes  as  the  imaginary  part  of  (u  + i v - u ) 1 , where  v ,,  is  the 

eff  D eft 

Ln  v1 

effective  collision  frequency,  u w — . rrj  , L is  the  density 

D * R V n 

e 

gradient  scale  length,  R is  the  major  radius  of  the  torus,  u*  is 
the  drift  frequency  and  is  the  electron  thermal  velocity  [8,  9,  12]. 


For  the  untrapped  particles,  the  real  part  of  F is  non-zero  only 

for  those  particles  with  v = u/k  . However,  for  drift  waves, 

u/k^  « , so  that  untrapped  runaway  electrons  are  not  affected  by 

the  magnetic  flutter.  For  the  trapped  electrons,  the  real  part  of  F 

can  result  either  from  u = or  else  from  the  ii/  ,,  term  in  the 

D ef  f 

resonant  denominator.  Since  the  wave  frequency  is  of  order  or 

less,  trapped  electrons  with  energy  larger  than  3 or  4 times  the 
thermal  energy  cannot  be  in  resonance.  Similarly,  as  electron  energy 
increases,  1/  decreases  so  that  the  real  part  of  F from  either 
resonance  or  collisions  is  very  small.  Thus  runaway  electrons  are  not 
resonant  with  any  wave  in  the  spectrum  and  are  not  preferentially 
diffused  by  the  magnetic  flutter  associated  with  drift  waves. 
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